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Abstract 

Dislocation based modeling of plasticity is one of the central challenges at the crossover of materials science 
and continuum mechanics. Developing a continuum theory of dislocations requires the solution of two long 
standing problems: {i) to represent dislocation kinematics in terms of a reasonable number of variables 
and {a) to derive averaged descriptions of the dislocation dynamics (i.e. material laws) in terms of these 
variables. The kinematic problem (j) was recently solved through the introduction of continuum dislocation 
dynamics (CDD), which provides kinematically consistent evolution equations of dislocation alignment ten¬ 
sors, presuming a given average dislocation velocity (Hochrainer (2015), Philos. Mag. 95 (12), 1321-1367). 
In the current paper we demonstrate how a free energy formulation may be used to solve the dynamic 
closure problem {ii) in CDD. We do so exemplarily for the lowest order CDD variant for curved dislocations 
in a single slip situation. In this case, a thermodynamically consistent average dislocation velocity is found 
to comprise five mesoscopic shear stress contributions. For a postulated free energy expression we identify 
among these stress contributions a back-stress term and a line-tension term, both of which have already been 
postulated for CDD. A new stress contribution occurs which is missing in earlier CDD models including the 
statistical continuum theory of straight parallel edge dislocations (Groma et al. (2003), Acta Mater. 51, 
1271-1281). Furthermore, two entirely new stress contributions arise from the curvature of dislocations. 
Keywords: Crystal plasticity. Continuum theory of dislocations. Thermodynamic consistency 


1. Introduction 

Plastic deformation of crystals is driven and controlled by the motion and interaction of dislocations. 
Dislocations are line like crystal defects which move when subjected to shear stresses and leave a permanent 
shear of the crystal when moving. Typically, metals contain a larger number of dislocations with line densities 
p in the range of 10^° to 10^^ m“^. These high numbers suggest that it should be possible to formulate a 
statistical mechanics theory of dislocations as a basis for crystal plasticity. However, as yet there is only 
a statistical theory available for strongly simplified systems of straight parallel edge dislocations (Groma 
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et al., 2003). This theory was originally developed from a Bogoliubov-Born-Green-Kirkwood-Yvon hierarchy 
adopted from interacting particle systems. For this derivation pair correlation functions have been obtained 
from statistical evaluation of a large number of quasi two-dimensional (2D) discrete dislocation simulations. 
Subsequently, part of the theory was recast into a phase-field type thermodynamic description (Groma 
et ah, 2007) derived from an effective free energy. A similar free energy for systems of straight parallel edge 
dislocations has recently been derived from a partition function by Kooiman et al. (2015). Transferring the 
statistical mechanics or thermodynamically inspired methods to systems of curved dislocations was doomed 
to failure as long as no suitable density variables capable of reflecting the kinematics of moving flexible lines 
were available. Such kinematically sound density variables have only recently been presented with the so 
called continuum dislocation dynamics (ODD) framework by Hochrainer et al. (2014) and Hochrainer (2015). 
A transfer of the statistical methods to curved dislocations seems impractical because of the difficulty to 
obtain sufficient statistics from three-dimensional discrete dislocation simulations. However, in the following 
we show that a thermodynamic approach is adoptable for curved dislocations. 

In the current paper we develop a thermodynamic formulation exemplarily for the lowest order GDD 
variant for curved dislocations in a single slip situation. The main result is that in this case a thermodynam¬ 
ically consistent average dislocation velocity is found to comprise five mesoscopic shear stress contributions 
as opposed to only one mesoscopic (back-) stress found in the quasi 2D case so far. For a postulated free 
energy expression we identify among these stress contributions a generalization of the quasi 2D back-stress 
term and a line-tension term, which have both been postulated in similar form for GDD before. Notably, 
a new stress contribution occurs which is missing in earlier GDD models including the quasi 2D case. Fur¬ 
thermore, two entirely new mesoscopic stress contributions are found to arise from the curved nature of 
dislocations. 

The paper is structured as follows: in Section 2 we briefly introduce notations and preliminaries from 
elasticity and crystal plasticity in small deformations. The quasi 2D theory of straight parallel edge dislo¬ 
cations serves both as point of departure and as reference in the sequel. We therefore recall the quasi 2D 
theory as derived from statistical mechanics and its thermodynamic reformulation in Section 3. In Section 
4 we briefly introduce the kinematics of the lowest order CDD for curved dislocations before we derive the 
general form of a thermodynamically consistent dislocation velocity v within this framework. Subsequently, 
we assume a free energy expression for the curved dislocation case in analogy to the one used in quasi 2D 
CDD. This free energy is used to illustrate the meaning of five mesoscopic stress contributions to the thermo¬ 
dynamically consistent dislocation velocity and to relate these stress contributions to expressions obtained 
in quasi 2D CDD and postulated for CDD of curved dislocations before. In Section 5 we summarize and 
discuss the results and depict open problems. 
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2. Notations and preliminaries 


We employ the convention that scalars are denoted with light face symbols while vectors and tensors are 
denoted with boldface symbols. Note, however, that we do not employ any notational difference between 
tensors of different orders, including vectors. Frequent use of coordinate notation shall avoid any ambiguities 
in this regard. We exclusively work in Cartesian coordinates with unit vectors i = 1, 2, 3}. Accordingly, 
we do not distinguish co- and contra-variant indices and employ the modified Einstein summation convention, 
where summation applies to all pairs of equal (lower) indices. The Kronecker delta (metric tensor) is denoted 
with and the totally antisymmetric Levi-Civita-symbol reads tijk- 

Partial derivatives in the coordinate directions shall be denoted with di\ the partial time derivative is 
usually denoted with 9*. We do not apply an according short hand notation for derivatives of functional 
kernels with respect to their arguments. 


2.1. Kinematics 

We restrict attention to small deformations such that the deformation is described by the displacement 
vector field u and the distortions by the displacement gradient Vm. In compliance with the small strain 
assumption we assume an additive split of the total distortion 

Vm = -k y3P‘, (2.1) 


into elastic and plastic distortion tensors, respectively. 

In crystal plasticity, plastic slip is confined to well defined slip planes which we characterize by slip plane 
normal n and slip direction m. The plastic distortion tensor in small deformations is additively composed 
of slip system specific shear tensors. In the sequel we only consider one active slip system. Therefore, the 
plastic distortion tensor is of the form 

ygP* = yn 0 m, (2.2) 


with 7 denoting the accumulated plastic slip on the slip system. We note that the slip direction m is the 
normalized Burgers vector b of dislocations on the slip system, i.e. m = b/b, where b denotes the modulus 
of the Burgers vector. 

The total strain tensor e derives from symmetrizing the displacement gradient 


~ 2 diUj) . 


(2.3) 


Accordingly, also the total strain is additively split into an elastic and a plastic strain, e = e®* + which 
each time arise from symmetrizing the respective distortion tensors. 
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2.2. Constitutive theory 


We assume linear elasticity such that the elastic constitutive equation connecting the stress tensor c 
with strain is 

(Tij = Ciju (dkUi - pfi'j = Cijkieti. (2.4) 

Volumetric forces will not be considered in this paper. Hence, conservation of momentum yields 


= 0 . 


(2.5) 


Due to conservation of moment of momentum the stress tensor is symmetric, aij = Gji. 

In the sequel we employ isotropic elasticity, such that we may write the elasticity tensor in terms of bulk 
modulus K and shear modulus G as 


^ijOkl 


Gijkl — H” G I ^ik^il d” ^il^ik 


ikOjl ~r OiiOjk 


JijOkl 


( 2 . 6 ) 


Note that the restriction to isotropic elasticity is convenient but not essential for the results of this paper. 


2.3. Boundary conditions and partial integration 

In the current paper we concentrate on bulk behavior and regard an infinite crystal which is loaded by 
a remotely applied external stress We frequently use partial integration, where we always assume that 
boundary terms vanish, e.g. by assuming that the involved variables vanish in the infinite or by employing 
periodic boundary conditions. In other words, we use 

JdJgdV = -Jfd,gdV (2.7) 

for suitable functions / and g. 


3. Continuum dislocation dynamics of straight parallel edge dislocations 

The statistical mechanics theory for systems of straight parallel edge dislocations developed by Groma 
et al. (2003) is as yet the only statistical continuum theory of dislocations derived from systematic averaging. 
The theory was successfully shown to produce essentially equivalent results to discrete dislocation simulations 
of straight edge dislocations treated as point particles with two possible signs in an elastic medium (Yefimov 
et ah, 2004). The theory henceforth serves as point of departure and as reference in the current paper. 

In the case of only edge dislocations we choose the coordinate system such that the Burgers vector points 
in 1-direction, b = 6ei, and positive edge dislocations point in 2-direction. The slip plane normal accordingly 
points in 3-direction, n = e^. Assuming additionally a plane strain state this dislocation system and its 
evolution can be completely described as quasi two-dimensional system in the 1-3 plane. 
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3.1. Continuum kinematics of straight parallel edge dislocations 

The dislocation system is described in terms of the so-called total dislocation density p and the net- 
dislocation density or geometrically necessary dislocation (GND) density k. The kinematics of plastic 
deformations as well as the evolution of the dislocation variables are described in terms of an average 
dislocation velocity v. For the regarded single slip system the plastic distortion tensor has the form = 
763 0 Cl, such that 7 = 0^^. The evolution of the plastic slip is obtained from Orowan’s equation. The 
evolution of the plastic variables consequently read (cf. Groma et al. (2003)) 

da = pvb, (3.1) 

dtp = -di {vk) , (3.2) 

dtK = -di {vp). (3.3) 

Before moving on to the constitutive equations we note that from the above one easily realizes the relation 
K = —dajh. We conclude that k defines the only non-vanishing component 0:21 = nh of the dislocation 
density tensor = EikidkPfj (Kroner, 1958). This justifies the denomination of k as the GND density. 

3.2. Constitutive theory for straight dislocations 
3.2.1. Results from statistical mechanics approach 

In order to employ the kinematic evolution equations (3.1)-(3.3) as a material law, a constitutive closure 
is needed which defines the average dislocation velocity v in terms of the current stress and dislocation state. 
The dislocation velocity was obtained from statistical averaging in the quasi 2D theory of straight parallel 
edge dislocations (Groma et ah, 2003; Yefimov et ah, 2004) in the form 

V = Mbsign (rnet) (iT-netl - Tf) , (3.4) 

where M is a dislocation mobility coefficient, the function sign returns the signature of its argument and 
(•) denote the Macaulay brackets, which return their argument if it is positive and zero else. The so-called 
flow-stress Tf defines a threshold (shear) stress such that the velocity is zero as long as the absolute value 
of the net-shear stress Tnet is smaller than the threshold. The flow stress is well known to obey the form 
Tf = aGb0f) (called the Taylor relation), with a dimensionless coefficient a typically found between 0.2 and 
0.5 (cf., e.g. Mecking and Kocks (1981)). This form for Tf could also be derived from rigorous averaging in 
the quasi 2D case (Groma et ah, 2003). The net shear stress was derived in the same work as 


0161 ■— 4 ” 06 0 )- (^■^) 

The appearing stress contributions are the external shear stress Tex = from remote boundary 

conditions independent of dislocations and the self-consistent stress Tsc, which accounts for the long range 
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stress field due to distributions of geometrically necessary dislocations. The self-consistent stress is obtained 
as 

Tsc(x):= J K{x')nnd{x,x')dV\ (3.6) 

where rind (^, denotes the shear stress induced at location x by a dislocation at position x' (see, e.g. Hirth 
and Lothe (1968)). Furthermore, there appears a mesoscopic stress contribution called the back-stress, which 
is of the form 

a a a to ' 7 ^ 

Tb =- OiK = -aidi7, (3.7) 

P P 

with a dimensionless constant B. Note that we use the sign convention for the back stress employed in 
Yefimov et al. (2004), which differs from the one in Groma et al. (2003). 

In the sequel we will abandon the distinction between the external stress and the self-consistent stress 
as they sum up to the macroscopic or mean field resolved shear tress r = Tex + t^c- The mean field shear 
stress is simply the resolved shear stress from the stress tensor tr which is obtained upon solving for the 
displacement field u in (2.4) subject to the solenoidality of the stress tensor, given the accumulated plastic 
slip 7 and remote boundary conditions cf., e.g. El-Azab (2000); Yefimov et al. (2004); Sandfeld et al. 
(2013). In other words, we find r = cr 2 i and rewrite (3.5) as 

Tnet:='r-rb. (3.8) 

We note that the use of the eigenstrain formulation instead of the convolution (3.6) for the self-consistent 
stress naturally allows for dropping the assumption of isotropic elasticity. 

3.2.2. Phase field type reformulation 

The net shear stress Tnet introduced above plays the role of a driving force for the evolution of the 
dislocation system. Originally derived from statistical averaging, it was later obtained in a phase field type 
description as a variational derivative of a free energy of the dislocation system. However, the flow stress Tf 
is not derivable in this way. It was therefore suggested by Groma et al. (2015) that the flow stress should be 
considered as a friction-type ingredient of a non-constant mobility function. With such a mobility function 
we assume the dislocation velocity to have the generic form 

V = M (Tnet, p)bTnef (3.9) 

In order for this form to be consistent with (3.4) the mobility function must be defined as 

M(r„et,p) = Mi^^^. (3.10) 

I '^net I 

Note that the mobility vanishes in the case of a vanishing net shear stress as it vanishes whenever |Tnet| < T[. 

The free energy ^ of the quasi 2D GDD system is composed of a ‘classical’ elastic energy 

tf'®* depending on boundary conditions and the ‘mean field’ elastic deformation e®* due to inhomogeneous 
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distributions of dislocations and a defect related density cf. Groma et al. (2015); Kooiman et al. (2015). 
The elastic energy S'®* will be discussed in Section 4.2 in the more general setting of curved dislocations which 
comprises the straight dislocation case as special case. For the defect related energy we adopt the form 
which was derived in Groma et al. (2007) for the quasi 2D case mostly from scaling arguments. With 
respect to the functional form, the free energy obtained by Kooiman et al. (2015) from direct averaging is 
not much different. We concentrate on the low GND approximation (|k| << p) for which the defect free 
energy according to Groma et al. (2007) assumes the form 




def 


'f'(p, Hi)dV 


Gh^ 




dK, 


(3.11) 


with dimensionless constants A and B and a reference dislocation density po- Groma et al. (2007) do not 
specify the reference density because the evolution of the system remains independent of it. 

The back-stress can be derived from the variational derivative of the defect free energy with respect to 
the net-dislocation density n as 


1 d'F{p,K) 


= BGbdi 




(3.12) 


where the latter approximation follows if additionally to |k| << p one assume that the gradients dip and din 
are of the same order of magnitude. As will be discussed in Section 4.2, the mean field shear stress r derives 
from the elastic part of the free energy. Consequently, the net shear stress Tnet = t — Tb may be obtained 
from the above mentioned free energy functional such that (3.9) may be interpreted as a mobility times a 
thermodynamics driving force as is common in non-equilibrium thermodynamics. However, the driving force 
appears to be incomplete in the sense that it does not involve a dependency on the variational derivative of 
the free energy density with respect to the total dislocation density p. That such a term should not appear 
has been argued for in Groma et al. (2007). But this was related to an argument put forth in Groma et al. 
(2003), where a similar term was rationalized away on the ground that dislocations of opposite character 
when subject to the same stress are bound to move with the same velocity in opposite directions. As we 
will see below, the additional term derived in the current paper does not rely on allowing for asymmetric 
dislocation velocities. It will turn out that the negligence of this term is not thermodynamically consistent 
unless the free energy density would either be independent of p or depending on p linearly. 


4. Continuum Dislocation Dynamics for curved dislocations 

Continuum dislocation dynamics comprises a whole hierarchy of kinematic theories obtainable by trun¬ 
cating a series of alignment tensors (Hochrainer, 2015) approximating a higher dimensional description 
(Hochrainer et al., 2007). In principle CDD is neither restricted to a single slip system nor to purely conser¬ 
vative dislocation motion. However, the only worked out examples so far are the two lowest order closures 
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for single slip systems subject to the assumption of an isotropic dislocation velocity (i.e., independent of 
dislocation character). In the sequel we will restrict attention to the lowest order truncation with internal 
variables: total dislocation density p, dislocation density vector p, and dislocation curvature density q. We 
note that the transfer of the thermodynamic approach presented in the following to higher order truncations 
and non-isotropic velocities is straight forward. 

4-.1. Kinematics of single slip CDD 

We employ the evolution equations in the conservative form for q as derived in Hochrainer (2015). The 
evolution of the plastic distortions and the dislocation variables then assume the form 


aft'' 

= pvbnimj, 

(4.1) 

dtp 

= di (vEijPj) + vq, 

(4.2) 

dtPi 

= -Sijdj (vp ), 

(4.3) 

dtq 

= di [vqi - Pjidjv) , 

(4.4) 


where = EikjUk denotes the operator performing a cross product with the slip plane normal n. Again, v 
denotes a dislocation velocity which is in the current framework assumed to be described by a scalar such 
that dislocation segments of different character will move with the same speed in directions perpendicular to 
their line-direction within the glide plane. The curvature vector qi and the second order dislocation density 
tensor py appearing in the evolution equation of q (4.4) are unknown and subject to closure assumptions. 
In the current work we do not specify the closure assumptions and refer to Monavari et al. (2014) and 
Hochrainer (2015) for discussions of this purely kinematic closure. We note that an impressive demonstration 
of the kinematic consistency of equations (4.1)-(4.4) was recently presented by Sandfeld and Po (2015) in a 
comparison with simplihed discrete dislocation simulations. 

The plastic distortions from a single slip system are again characterized by the accumulated plastic slip 
7 , such that = yn (8) m. The role of the GND density k in the quasi 2D theory is now adopted by the 
dislocation density (or GND) vector which relates to the accumulated plastic slip through 

7 1 

Pi — E,ijf^dj—Tlk — ’ (^■^) 

We note that the kinematic evolution equations for straight parallel edge dislocations (3.1)-(3.3) arise as 
special case of the above evolution equations (4.1)-(4.4) if all quantities are constant in edge direction and 
the curvature density q vanishes identically. With the coordinate system introduced in Section 3.2 the GND 
density k is in this case the 2-component (i.e. edge component) of the dislocation density vector, k = p 2 - 

4-2. Thermodynamically consistent single slip CDD 

While there is currently no free energy available for GDD of curved dislocations it seems natural that a 
free energy •P of the CDD system will also be composed of an elastic energy depending on the elastic 


strain e®* due to boundary conditions and the ‘mean field’ deformation due to inhomogeneous distributions of 
dislocations, and a defect related energy Both energy contributions arise from integrals over according 
energy densities which we define through respective terms in the following equation 

= I : c : e^'d^ + f p , q)dV. (4.6) 


In the sequel we derive a condition for the functional form of the dislocation velocity as a function 
of the thermodynamic forces defined through the variational derivatives of the free-energy functional. The 
development of CDD mainly targets at face centered cubic (fee) crystals, where -aside from high temperature 
applications- thermal effects are usually assumed to be negligible (Groma et ah, 2007). We therefore also 
ignore temperature effects and consequently find the requirement of thermodynamic consistency such that 
the spontaneous evolution of the system may not increase the free energy. In other words, we require that 


dtf' 

— < 
dt 


0 . 


(4.7) 


Assuming that the free energy does not depend explicitly on time the total time derivative is obtained from 
the chain rule as 

d<f- _ /• 

dt J 

For the elastic part we find using partial integration (2.7) and (4.1) 


dp dt 


dll' dpi 
dpi dt 


d'l' dq 
dq dt 


dy. 


(4.8) 


dt 


cle- 




J (didtUj - (JijdV 

J —dtUjdiaij — pvbaijUimjdV 

/ —pvbrdV. 


(4.9) 


The last identity follows from the solenoidality of the stress tensor (2.5) and we introduced the resolved 
(mean-field) shear stress r = aijUirrij. 

For the explicitly defect related part (note that the mean field stress depends on the GND density via 
the incompatibility) we find from (4.8) and repeated application of partial integration 


d«? 


def 


dt 


[d, {vEijPj) + ^^9] - {vp)] + ^ [di {vqi - Pjidjv)] I dl/ 


d'P 
dp 

^d'l' d'^ ( dW\ d'P ( d'P 

Qp+^Qp+ P^3 d, i^PjA 


dF 


Inserting (4.10) and (4.9) into (4.8) the requirement for thermodynamic consistency reads 

f.d'F d'F f d<I'\ d'F f ^ x ..r ^ r, 

phr SijPjOi Qp ^ ^ Qp ^ Qp] Qq ^3 I 


(4.10) 


(4.11) 
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A thermodynamically consistent dislocation velocity must therefore have the opposite sign than the term in 
square brackets. This is ensured if the velocity contains the negative multiple of the term in brackets as a 
factor. We assume the velocity has the same form as in the quasi 2D case, 

V = M{Tmu p)bTnet- (4-12) 


Moreover, we may safely assume that the flow stress in the case of curved dislocations also obeys the Taylor 
relation Tf = aGh^fp such that the nonlinear mobility function M (rnet, p) may again be defined by (3.10). 

With the velocity of the form (4.12), thermodynamic consistency is ensured if the net shear stress Tnet is 
a negative multiple of the term in brackets in (4.11). In order to obtain terms which have the dimension of 
stress we divide the term in brackets by pb and arrive at 


'^net 


pb "dp 


q d'l' 
pb dp 




pb " dq 


7b^^ 



(4.13) 


4.2.1. Thermodynamic consistency from a postulated energy functional 

In order to understand the appearing mesoscopic stress contributions in (4.13) and to relate them to 
the terms obtained in the straight dislocation case we now assume the free energy density in the curved 
dislocation case in analogous form to the one for straight dislocations as 


'^{P.TP) = Gb"^ Ap\n 


1 B.jpipj ^ ^ 1 

2 p ^ 2 


(4.14) 


with dimensionless constants A, Bij and C and a reference density po • In a spatially and elastically isotropic 
theory as developed here the tensor B is expected to be spherical, i.e Bij = B5ij, with a constant B. 
Assuming this form we find 


^ = GV 


dT> 

~^P 

dT 

dpi 

dT> 

dq 


A f In ( 1 + 1 1 - 

Po 


B piPi 

2 p2 


PM 


= Gb'^B— 
P 

= Gb^G\. 


(4.15) 

(4.16) 

(4.17) 


Before we proceed with inserting these expressions into the definition of the net shear stress (4.13), we 
make some simplifying assumptions. Please note that these approximations are mostly done here to arrive 
at readable and partly familiar expressions. In principle any of the following negligences may violate ther¬ 
modynamic consistency in special situations. As in the quasi 2D theory we assume small GND densities, 
i/PiPi << p, and additionally we assume that the average radius of curvature r = p/q is much larger than 
the average dislocation spacing, i.e. r >> If^/p. Because the curvature density is g = pjr this implies 
« p^, such that we approximate (4.15) by 


In ( ^ ) + 1 
Po 


dT 


Gb'^A 
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(4.18) 







With regard to the spatial derivatives of the variational derivatives in (4.13) we again assume that the 
spatial derivatives of all density variables p,p, and q are of the same order of magnitude. We then obtain 
the approximate expressions 



a AGb^^, 

dp 

P 

di' 

e BGb^^^^^ 

dpj 

P 

d'B 

a CGb^^. 

* dq 



Inserting these approximate expression into the net shear stress (4.13) we obtain 


Tnet ^T + Gb 


- Y^d^p - 

P P 


In 


+ 1 - 


Be 


Cq, 


— djPi + ^^^diq ' 
p p-^ 


c 

p 


5, 


PP 


diq 


(4.19) 

(4.20) 

(4.21) 

(4.22) 


We interpret each of the summands in the curly braces (upon multiplication with Gb) as a mesoscopic shear 
stress contribution. We begin with the term which generalizes the back stress found in the theory of straight 
parallel edge dislocations, which is the third term in the braces. That this is the back stress seems obvious 
because the term involves a gradient of the GND density. Because the GND density is itself a gradient of 
the plastic slip, this term is again a second order strain gradient and we find 


BGb, 




Th = 


BG 

P 


(4.23) 


The last expression in (4.23) contains the Laplace operator within the glide plane (5^- — nitij) didj. The 
according identity in (4.23) is obtained from the antisymmetry Sji = —Sij and the fact that £ik£kj = 
— (5y — niUj). The form of the back stress term is remarkable not only because it naturally generalizes 
the expression obtained in the quasi 2D theory, but because the Laplace operator of the plastic shear has 
been postulated in phenomenological strain gradient models used to account for strain gradient effects in 
small scale plasticity (Aifantis, 1987). Notably, however, the length scale in front of the current back stress 
expression is consistently found as the only characteristic length scale in a dislocation system, which is the 
average dislocation spacing governed by the dislocation density. 

A further remarkable term is the first term in the curly braces in (4.22), which we call the (density) 
gradient stress 


AGbSijPj ^ ^ 
Tg — ^2 b)ip. 


(4.24) 


This shear stress contribution is of special interest for two reasons: (?) it contains the gradient of the total 
dislocation density and (??) the term stems from the flux part of the total dislocation density evolution and 
should therefore also appear in a 2D theory of straight parallel edge dislocations. For realizing the latter 
point we note that the term is a contraction of the density vector p tilted by 90° in the glide plane and the 
gradient of p. This term therefore reads 


AGh „ 

-5-Kdip 
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(4.25) 









in the straight dislocation case, where we remind that k = p 2 - That the term contains a gradient of the 
total dislocation density at first sight seems unphysical, as it appears to drive a diffusive behavior for the 
total dislocation density. However, the gradient stress additionally requires the existence of GNDs, which 
on the one hand suggests the appearance of a local stress field and on the other hand renders the expression 
directional such that the contribution is not actually diffusional. 

The remaining mesoscopic stress contributions are unique to the theory of curved dislocations. One of 
them has been postulated earlier by Hochrainer et al. (2014) as a line tension term. To see the connection 
we introduce the average dislocation curvature (the inverse of the average radius r discussed before) which 
is defined as k = q/p. The line tension term has been postulated to be of the form rit = Tk with a line 
tension T oc Gb. The according term in (4.22) is the second one in the curly braces which reads 


rit = AGh (^In + 1^ k. (4.26) 

The current form of the term is irritating because it explicitly depends on the reference density po which is 
not supposed to appear in the evolution equations according to Groma et al. (2007). We shall postpone the 
discussion of this point to Section 5, and only note at this point, that regardless of the specific form of the 
free energy, a line-tension type shear stress will appear in thermodynamically consistent constitutive laws 
unless the free energy would be independent of p. 

The two stress contributions which contain derivatives of the curvature density, 


Tql = 


GGbq^ 


diq, 


and 


GGb 



(4.27) 


are completely new and unique to the theory of curved dislocations. As yet we cannot offer physical 
interpretation of these terms. But we observe that Tqi is structurally similar to Tg. Note that Tqi is a 
contraction of the curvature vector q and the gradient of the scalar curvature density q similar to Tg being a 
contraction of the dislocation density vector p and the gradient of the density p. The second stress term can 
be further transformed by employing the fact that the dislocation curvature vector arises as the divergence 
of the second order alignment tensor, i.e qi = djPji (Hochrainer, 2015). Accordingly, we obtain 


Tq2 — Tql + 


CGbp,. ^ 



(4.28) 


Because the last term loosely resembles the back stress (in being a second order derivative) we denote it as 


— 





(4.29) 


For harmonizing the appearance we define the gradient curvature stress as Tqg := 2Tqi and finally obtain the 
net shear stress as 

Tnet = T -l-Tg - Tit - Tb -I- Tqg -l-Tqb. (4.30) 
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We derived the interpretations of the stress contributions for the free energy density (4.14) postulated in 
analogy to the free energy expression presented by Groma et al. (2007) in the quasi 2D case. But we may 
transfer these interpretations to any free-energy density, such that we define in general 


Tit := 


pb "dp' 

kd'P 
b dp ’ 


1 


Tb 


d'P 


'qg 


'^qb • — 


2«i „ a® 


Pji 


di d, 


d^ 


(4.31) 

(4.32) 

(4.33) 

(4.34) 

(4.35) 


pb ^ y'" dq , 

With these definitions entered into the net shear stress (4.30) the velocity may also in general be defined by 
(4.12). 

We finally find that a thermodynamic consistent dislocation velocity of the form (3.9) for the quasi 2D 
case requires the net shear stress to be of the form 


2D 


T “T Tg 5 


(4.36) 


with 


^2D 


^2D 


1 dw 


(4.37) 

(4.38) 


Taking the quasi 2D free energy from (3.11) the above definitions, in conjunction with the simplifying 
assumptions employed throughout the current Section, reproduce (3.12) for and (4.25) for r|°. Notably, 
we conclude that the statistical mechanics theory developed by Groma et al. (2003) is not thermodynamically 
consistent with the free energy proposed in Groma et al. (2007). But for the small GND case treated in the 
named papers we would expect that the gradient stress (4.25) would be much smaller than the back stress 
(3.12), at least when the constants A and B in the free energy (3.11) are of the same order of magnitude. 


5. Discussion 

In the current paper we present the first thermodynamically consistent phase field-like framework for 
GDD comprising the variants for straight parallel edge dislocations and the single slip theory based on the 
internal state variables: total dislocation density p, dislocation density vector p, and dislocation curvature 
density q. The constitutive equations for the dislocation velocity in ODD had earlier been assumed in 
analogy to the quasi 2D case which made them appear somewhat ad hoc because they were neither obtained 
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from statistical averaging nor from a general principle in the spirit of non-equilibrium thermodynamics. In 
the current contribution we showed how this deficiency may be overcome, once a free energy is formulated 
in terms of the internal state variables of CDD. The results obtained for a free energy functional postulated 
in analogy to a free energy functional available for straight parallel edge dislocations are quite satisfactory 
because they justify the form of a back stress term and the appearance of a line tension term, both of 
which have been postulated in similar form for CDD before (Zaiser et ah, 2007; Hochrainer et ah, 2014). 
Surprisingly, thermodynamic consistency requires the appearance of a so-called gradient stress which is 
missing in the theory of Groma and co-workers both in the statistical mechanics as in the phase field 
approach to quasi 2D dislocation systems. With the current derivation we not only find that this term is 
missing in the straight dislocation case but also derived the vectorial generalization of its definition for CDD 
of curved dislocations. Besides the line tension term, two additional shear stress contributions were obtained 
which relate to the curvature of dislocations. These terms are as yet not well understood and it seems too 
early to analyze them in detail before more is known about the actual dependence of the free energy on the 
curvature density q. 

The line tension term th in (4.26) was found to be depending on the reference density po used in the 
logarithmic expression in the energy density. While Groma et al. (2007) would consider this as contradicting 
a scale invariant theory, we note that an energy contribution of the form pin (p/po) has been derived earlier, 
e.g. by Wilkens (1969) and Berdichevsky (2006) which both assign a physical meaning to the reference 
density po. In Wilkens (1969) the energy of a system of dislocations with density p is reported to have the 
form pin (i?e/?'o), with an outer cut-off radius i?e and an inner cut-off radius rp, in analogy to the energy 
of single dislocations (Birth and Lothe, 1968). For certain random distributions of dislocations the outer 
cut-off radius is found to scale like the average dislocation spacing, i?e oc l/i/p, due to screening effects. 
The inner cut-off radius is taken as the dislocation core radius which is on the order of the magnitude of the 
Burgers vector rp ~ b. We may therefore transform the energy expression of Wilkens (1969) such that 

pln(i?e/ro) = Bpln(p/5“^) , (5.1) 

with a (negative) constant A, which contains all unknown factors. The reference density in this case is 
Po = 6“^, which is well defined and an obvious upper bound or saturation value for admissible dislocation 
densities. That pp should be an upper bound for the dislocation density was also put forward by Berdichevsky 
(2006) who suggested the same functional form for the energetic contribution of the total dislocation density. 
Berdichevsky (2006) interprets pp as a saturation density of the crystal, and obtains for Al and Ni the value 
Po ~ 3 • 10^“^ m“^. We shall not dwell upon the question of the correct definition of po at this point. The 
important message to take away is that the reference density is not an arbitrary scaling parameter which 
has to drop out of the evolution equations, but a well defined physical quantity. As was pointed out to me 
by I. Groma, the quasi 2D system is scale free due to the 1/r-dependence of the dislocation interactions, 
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while in 3D systems dislocation core effects as the line tension may well introduce a further parameter with 
length dimension, as e.g., the core radius. This line of thought is in accordance with the above cited work 
by Wilkens. 

As the combination of flow stress and net shear stress in CDD essentially translates into the flow rule 
for phenomenological plasticity laws it is also remarkable that the back stress expression obtained in the 
current paper resembles a second order strain gradient expression postulated by Aifantis (1987) in a phe¬ 
nomenological theory before. An important difference to the latter theory is, however, that the prefactor, 
which is required to have the dimension of a length squared, is not a constant but 1/p, that is the square of 
the average dislocation spacing l/y^, which is the only characteristic length scale of the system. The newly 
found gradient stress additionally yields a dependence of the ‘flow rule’ on the first order strain gradient, i.e. 
the GND tensor. While such dependencies have been postulated numerous times in strain gradient plasticity 
since the pioneering works of Ashby (1970); Fleck and Hutchinson (1993); Nix and Gao (1998); Gao et al. 
(1999) to name just a few, GNDs were usually assumed to increase the flow stress. By contrast, the gradient 
stress derived in the current work is a signed quantity, which may obstruct or facilitate plastic deformation 
depending on the relative direction of the gradient of plastic slip and the gradient of total dislocation density. 
Furthermore, the required coefficient is again not a constant but involves the total dislocation density and 
its gradient. 

The herein derived theory of course remains preliminary in many respects. Most obviously, the free 
energy functional is not yet known for the considered single slip GDD theory. But this problem is common 
to all available strain gradient plasticity theories, as well. Note, however, that GDD is solely based on 
physically defined internal variables such that there is reason to expect that an according free energy density 
is derivable from dislocation theory . In fact, during the review process Michael Zaiser proposed a framework 
for deriving the free energy for systems of curved dislocations in terms of the dislocation density vector and 
the second order dislocation alignment tensor using pair correlation functions (Zaiser, 2015). Another open 
problem is the form of the mobility function which remains on the level of an educated guess so far. As 
for the flow stress we note that the Taylor relation is a universal finding at least in fee crystals. Because 
the dimensionless prefactor a is known to fall into a narrow range, 0.2 < a < 0.5 and to be only weakly 
depending on details of the dislocation structure it seems feasible to treat it as a htting parameter rather 
than trying to obtain it from ‘first principles’. Of course, the mobility may also be defined to yield a smoother 
transition to plastic flow around the flow stress. 

Another point obviously missing in the current paper is the consideration of boundary conditions. The 
decision to neglect boundary conditions was mostly done to keep the paper focused. The derivation of 
boundary conditions is possible through the application of Gauss’ integration theorem in place of the partial 
integration (2.7) in the time evolution of the free energy. 

In the current derivation we solely developed a single slip theory. But of course any crystal plasticity 
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theory needs to be designed for multiple active slip systems. The most simple option would be to treat 
multiple slip mostly as a superposition of single slips ~ with the exception that the flow stress needs to 
consider dislocations on other slip systems s, e.g. in the generalized Taylor relation (Franciosi et ah, 1980) 
T® = ass'Ps', with dimensionless parameters ass' characterizing the ‘frictional interaction strength’ 

between slip systems s and s'. Truly multiple slip theories will of course need to consider shear stress 
contributions from correlations between dislocations on different slip systems. The derivation of such a 
theory will require further research on the kinematic as well as on the energetic/kinetic side. 

In the present paper we concentrated on the lowest order CDD theory for curved dislocations. In 
Hochrainer (2015) it was demonstrated how more refined single slip CDD theories may be systematically 
derived from an alignment tensor expansion of a higher dimensional dislocation theory. In the very paper 
two extensions are sketched, which is (i) a theory again employing an isotropic dislocation velocity v but 
considering the evolution of the second order dislocation alignment tensor pij , and (ii) a theory for anisotropic 
velocities described by a second order velocity tensor Vij. Assuming that the dependency of the free energy 
on higher order tensors is known, as suggested by the recent work of Zaiser (2015), the derivation of 
a thermodynamically consistent dislocation velocity put forward in this paper can be straight-forwardly 
applied to these more refined theories. Maybe the most surprising finding of the current paper is that the 
admissible and always size-dependent kinetic equations are essentially dictated by the kinematic equations 
of CDD. 
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